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In this work we consider the steady state scaUng behavior of directed 
percolation around the upper critical dimension. In particular we 
determine numerically the order parameter, its fluctuations as well 
as the susceptibility as a function of the control parameter and the 
conjugated field. Additionally to the universal scaling functions, 
several universal amplitude combinations are considered. We 
compare our results with those of a renormalization group approach. 
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1 INTRODUCTION 

The concept of universality is one of the most impressive features of continuous 
phase transitions. It allows to group the great variety of models into a small 
number of universality classes (see for a recent review). Models within one 
class share the same critical exponents. Furthermore their corresponding scaling 
functions become identical close to the critical point. Often, universality classes are 
also characterized by certain amplitude combinations, which are merely particular 
values of the scaling functions. The most prominent examples of universal behavior 
are the coexistence curve of liquid-vapor systems and the equation of state in 
ferromagnetic systems (e.g. ^ Deciding on a systems universality class by 
considering the scaling functions instead of critical exponents appears to be less 
prone to errors in most cases. While for the latter ones the variations between 
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different universal classes are often small, the amplitude combinations and therefore 
the scaling functions may differ significantly (see |^). 

Wilsons renormalization group (RG) approach lays the foundation for an under- 
standing of universality It also yields a tool for computing critical exponents 
as well as the universal scaling functions. While critical exponents emerge from 
local properties near a given fixed point, scaling functions require the knowledge 
of the full RG flow along the trajectories between neighboring fixed points. This 
illustrates why scaling functions are more sensitive than the corresponding expo- 
nents. 

The RG explains the existence of an upper critical dimension above which 
the mean-field theory applies, i.e., classical theories, which neglect strong fluctua- 
tions and correlations, provide correct estimates to the critical exponents and scal- 
ing functions. Below fluctuations become relevant and the mean-field scenario 
breaks down. At the upper critical dimension the RG equations yield mean-field 
exponents with logarithmic corrections [3 EI • 

In contrast to equilibrium critical phenomena less is known in the case of non- 
equilibrium phase transitions. This is due to the fact that a generalized treatment 
is not possible, lacking an analog to the equilibrium free energy. The rich and of- 
ten surprising variety of phenomena has to be studied for each system individually. 
The scaling behavior of directed percolation (DP) is recognized as the paradigm of 
the critical behavior of several non-equilibrium systems which exhibit a continuous 
phase transition from an active to an absorbing state (see e.g. [S]). The widespread 
occurrence of such models describing phenomena in physics, biology, as well as cat- 
alytic chemical reactions is refiected by the well known universality hypothesis of 
Janssen and Grassberger: Short-range interacting models, which exhibit a contin- 
uous phase transition into a unique absorbing state belong to the DP universality 
class, provided they have a one-component order parameter and no additional sym- 
metries |in(lllj. Different universality classes are expected to occur in the presence 
of additional symmetries, like particle conservation J2j, particle- hole symmetry 
(compact directed percolation) ,13,, or parity conservation (e.g. branching annihi- 
lating random walks with an even number of offsprings Jl]). Other model details, 
such as e.g. the geometry or shape of a lattice, are expected to have no influence 
on the scaling behavior in the vicinity of the critical point. 

The universality hypothesis still awaits a rigorous proof. Amazingly, numerous 
simulations suggest that the DP universality class is even larger than expected. 
It turns out that the hypothesis defines only a sufficient condition but fails to 
describe the DP class in full generality (see 0| for a detailed discussion). For in- 
stance, the pair contact process (PGP) is one of the simplest models with infinitely 
many absorbing states exhibiting a continuous phase transition jJS]. It was shown 
that the critical scaling behavior of the one-dimensional PGP is characterized by 
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the same critical exponents |15( I16j as well as by the same universal scaling func- 
tions as DP Thus despite the different structure of the absorbing phase the 
one-dimensional PCP belongs to the DP universality class. This numerical evi- 
dence confirms a corresponding RG-conjecture |TS]. But one has to mention that 
a recently performed RG analysis conjectures a different scaling behavior of both 
models in higher dimensions jl9j . 

In this work we consider the universal scaling behavior of directed percolation 
in various dimensions. Whereas most investigations on DP follow the seminal work 
ref. 120] and thus focus on activity spreading we examine the steady state scaling 
behavior for D > 2. We determine the universal scaling functions of the order 
parameter (i.e. the equation of state) and its fluctuations. Furthermore we consider 
certain universal amplitude combinations which are related to the order parameter 
and its susceptibility. These amplitude combinations are immediately related to 
particular values of the universal scaling functions and are of great experimental 
interest j^. We will see that the numerically obtained universal scaling functions 
and the related universal amplitude combinations allow a quantitative test of RG- 
results. The powerful and versatile e-expansion provide estimates of almost all 
quantities of interest, e.g. the critical exponents and the scaling functions (see 
e.g. |21j). Unfortunately it is impossible to estimate within this approximation 
scheme the corresponding error-bars. Thus it is intriguing to compare our results 
with those of RG analysis [121 HB] • 

Furthermore we focus on the phase transition at the upper critical dimension 
= 4. There the usual power-laws are modified by logarithmic corrections. 
These logarithmic corrections are well established in equilibrium critical phenom- 
ena [71 IH] but they have been largely ignored for non-equilibrium phase transitions. 
Due to the considerable numerical effort, sufficiently accurate simulation data for 
non-equilibrium systems became available only recently: Investigated systems in- 
clude self-avoiding random walks [SlISSI, self-organized critical systems |26l l27]. 
depinning-transitions in disordered media [SH], isotropic percolation [22], as well 
as absorbing phase transitions [32] • On the other hand, the numerical advance 
triggered further analytical RG calculations yielding estimates for the logarithmic 
correction exponents for the respective systems [HTl 1^ 15^ OH] . 

The outline of the present paper is as follows: The next section contains the 
model definition and a description of the method of numerical analysis. In Sec. [HI 
we describe the scaling behavior at the critical point and introduce the critical 
exponents as well as the universal scaling functions. The numerical data of the order 
parameter and its fluctuations are analyzed in Sec.[l]below (D = 2, 3), above {D = 
5), and at the upper critical dimension (D = 4). Several amplitude combinations 
are considered in Sec. [3 Concluding remarks are given in Sec. [HI 
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2 MODEL AND SIMULATIONS 

In order to examine the scaling behavior of the D-dimensional DP universaUty 
class we consider the directed site percolation process using a generalized Domany- 
Kinzel automaton ^Ij. It is defined on a D + 1-dimensional body centered cubic 
(bcc) lattice (where time corresponds to the [0,0,..., 0, 1] direction) and evolves 
by parallel update according to the following rules: A site at time t is occupied 
with probability p if at least one of its 2^ backward neighboring sites (time t — 1) 
is occupied. Otherwise the site remains empty. Furthermore, spontaneous particle 
creation may take place at all sites with probability Pq. Directed site percolation 
corresponds to the choice po = 0. The propagation probability p is the control 
parameter of the phase transition, i.e., below a critical value p^ the activity ceases 
and the system is trapped forever in the absorbing state (empty lattice). On the 
other hand a non-zero density of (active) particles is found for p > p^. The 
best estimates of the critical value of directed site percolation on bcc lattices are 
p, = 0.705489(4) ^ for: D = 1 and p, = 0.34457(1) W for D = 2. 

The order parameter of the absorbing phase transition vanishes at the critical 
point according to 

Pa OC (1) 

with 5p = {p — Pc)/Pc- Furthermore the order parameter fluctuations Ap^ = 




Figure 1: The two-dimensional directed percolation order parameter p^ as a func- 
tion of the particle density for zero field (symbols) and for various values of the 
external field {h = 310'^, 10"^, 210'^, 510"^, 10"^) (lines). The inset displays 
the order parameter fluctuations Ap^ for zero field (symbols) and for various values 
of the external field h (lines). 
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L^HpD - iP^f) diverge as 

Ap^ oc 5p"^'. (2) 

The fluctuation exponent 7' obeys the scahng relation 7' = Di'j^ — 2f3 16 where i^x 
describes the divergence of the spatial correlation length at the critical point. The 
critical behavior of the order parameter is shown in Fig. ^ for D = 2. The data are 
obtained from numerical simulations of systems with periodic boundary conditions. 
Considering various system sizes L we take care that our results are not affected by 
finite-size effects. The system is started from a random initial configuration. After 
a certain transient regime a steady state is reached, which is characterized by the 
average particle density and its fluctuations Ap^. 

Similar to equilibrium phase transitions it is possible in DP to apply an external 
field h that is conjugated to the order parameter. Being a conjugated field its has 
to destroy the absorbing phase and the corresponding linear response function has 
to diverge at the critical point i.e., 

Xa = ^ - 00. (3) 

In DP the external field is implemented jHlEl a-s a spontaneous creation of particles 
(i.e. po = h > 0). Clearly, the absorbing state and thus the phase transition 
are destroyed. Figure^ shows how the external field results in a smoothening of 
the zero-field order parameter curve. The inset displays that the fluctuations are 
peaked for finite fields. Approaching the transition point (h 0) this peak becomes 
a divergence signalling the critical point. 

3 UNIVERSAL SCALING FORMS 

Sufficiently close to the critical point the order parameter, its fluctuations, as well 
as the order parameter susceptibility can be described by generalized homogeneous 
functions 

p^i6p,h) ~ A-^ ii(a/pA,a,/iA"), (4) 
a^Ap^{6p,h) ~ X^' D{ap6p X,ahhX'') , (5) 

a^Xa(^P7^) ~ C{ap6p X,a,^hX'^) . (6) 

Note that these scaling forms are valid for D ^ D^. At the upper critical dimension 
they have to be modifled by logarithmic corrections j3Uj . Taking into considera- 
tion that the susceptibility is deflned as the derivative of the order parameter with 
respect to the conjugated fleld [Eq. we flnd C{x,y) = dyR{x,y), = a^^, as 
well as the Widom scaling law 

j = a-p. (7) 
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The universal scaling functions R, D, and C are identical for all models be- 
longing to a given universality class whereas all non-universal system-dependent 
details (e.g. the lattice structure, range of interactions, the update scheme, etc.) 
are contained in the so-called non- universal metric factors a^, a^, and Ca The 
universal scaling functions can be normalized by the conditions 0) = -R(0, 1) = 
D{0, 1) = 1. In that case the non-universal metric factors are determined by the 
amplitudes of the corresponding power-laws 

p^i6p,h = 0) ~ {a,6pf, (8) 
p^{6p = 0,h) ~ {a,hf/% (9) 
a^Ap^{6p = 0,h) ~ (a,/i)-^>. (10) 

Furthermore we just mention that C(0, 1) = f3/a, following trivially from the defi- 
nition of the susceptibility. 

Usually, an analytical expression for the scaling functions is only known above 
D^, where mean-field theories apply. In the case of directed percolation the mean- 
field scaling functions are given by (see e.g. [SH|) 



RMF{x,y) 
Cm¥{xi y) 

i.e., the mean-field exponents are /3mf = l, o"mf = 2, 7mf = 1, and = 
(corresponding to a finite jump of the fiuctuations). Below the universal scaling 
functions depend on dimensionality and are unknown due to a lack of analytical 

Table 1: The non- universal quantities for site directed percolation on a bcc lattice 
for various dimensions. The uncertainty of the metric factors is less than 7%. The 
values for D = 1 are obtained from a previous work [17 . 



D 


Pc 




ah 


Ca 


1 


0.705489 ± 0.000004 


2.498 


0.114 


9.382 


2 


0.344575 ± 0.000015 


0.795 


0.186 


9.016 


3 


0.160950 ±0.000030 


0.417 


0.328 


11.91 


4 


0.075582 ± 0.000017 


14.70 3.055 


59.80 


19.19 


5 


0.035967 ± 0.000023 


0.114 


0.174 


42.49 



/ / \ 2 

X / ( X\ 



2 + 



\ly + {x/2f 
1 



2 + (x/2) 



(11) 

(12) 
(13) 
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Table 2: The critical exponents of directed percolation for various dimensions D. 
The one-dimensional values were obtained in a famous series expansion by 
Jensen ^39 . For D = 2 and D = 3 the authors investigated activity spreading 
and the presented exponents are derived via scaling relations. A complete list of 
all critical exponents of DP can be found in 9 . The symbol * denotes logarithmic 
corrections to the power-law behavior. 

D IW 2110] 301] 4MF 

P 0.276486(8) 0.584(4) 0.81(1) 1* 1 
a 2.554216(13) 2.18(1) 2.04(2) 2* 2 
7' 0.543882(16) 0.300(11) 0.123(25) 0* 

solutions. In this case the scaling functions have to be determined numerically or 
via approximation schemes, e.g. series expansions or e-expansion of RG approaches. 

In case of the mean-field solution {'y'^p = 0) the scaling form of the fluctuations 
[Eq. Q] reduces to 

a^Ap^{6p,h) ~ D{aj,5p X^a^hX"). (14) 

Some interesting properties of the universal scaling function D can be derived from 
this form. The non-universal metric factor Oa is determined by 

= A/>a(5p = 0,/.) 

using that D(0, 1) = 1. The value I?(1,0) is then given by 

Finally, it is worth mentioning that the mean-field scaling function D fulfills the 
symmetries 

D{l,x) = Dix-^/",!) (17) 

D{x,l) = D{l,x''') (18) 

for all positive x. In particular we obtain for x — > D{1,0) = D{oo,l) and 
l)(0, 1) = l)(l,oo), respectively. 

4 EQUATION OF STATE AND FLUCTUATIONS 

4.1 Below the upper critical dimension 

The scaling forms Eqs. H4I6|) imply that curves corresponding to different values of 
the conjugated field collapse to the universal functions R{x,l), D{x,l), C{x,l), 
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Figure 2: The universal scaling plots of the order parameter and its fluctuations 
(inset) for D = 2. The dashed line corresponds to an e-expansion of a RG ap- 
proach [221 ■ 

if p^{ahh)~^/^ , a£^/S.p^{a,^h)'^' , and a^Xs^io^hh)^^^ are considered as functions of 
the rescaled control parameter a^f^p (0^/1)"^/°". In a first step, the non-universal 
metric factors a^, a^, Oa are obtained from measuring the power-laws Eqs. H8I1U|) 
(see Table^). Here, the best known estimates for critical exponents, as given in 
TableI21 are used. 

Subsequently, the rescaled order parameter and its fluctuations as a function 
of the rescaled control parameter are plotted for two- and three-dimensional DP 
(Figs. 12121). A convincing data collapse is achieved, confirming the scaling ansatz 
as well as the values of the critical exponents. 

Besides the universal scaling function R(x, 1) the corresponding curve of an 
e-expansion obtained from a renormalization group analysis is shown in Figs. 12131 
Using the parametric representation j42l I43j of the absorbing phase transition, 
Janssen et al. showed that the equation of state is given by the remarkably simple 
scaling function j22j 

H{e) = e{2-e) + o{e^), (19) 

where e denotes the distance to the upper critical dimension = 4, i.e., e = D^—D. 
Here the scaling behavior of the quantities ^P^ and h is transformed to the 
variables R and 9 through the relations 



h5p = R(l-e), p. = R^-. 



(20) 
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Figure 3: The universal scaling plots of the order parameter and its fluctuations 
(inset) for D = 3 and for various values of the external field {h = 310~^, 10~^, 
210~^, 410~®, 510^^). The dashed line corresponds to an e-expansion of a RG 
approach 

The equation of state is given by 



with the metric factors a and b. The whole phase diagram is described by the 
parameter range R > and 9 G [0,2]. In Fig. l2l3l a comparison between the 
numerically obtained scaling functions and the analytical result of Eqs. ()19I21|) 
is made. The RG-data differ slightly from the universal function. As expected 
the differences decrease with increasing dimension and are especially strong in 
D = 1 This point is discussed in detail below. 

4.2 Above the upper critical dimension 

Above the upper critical dimension the scaling behavior of a phase transition equals 
the scaling behavior of the corresponding mean-field solution [Eqs. H11II13() ]. Plot- 
ting Pa/\/oh^ as a function of a^dp/ yja^h, the numerical data should collapse to 
the universal scaling function 



with the scaling variable x = a^bp j \J h. In Fig.^ we plot the corresponding 
rescaled data of the five-dimensional model. A perfect collapse of the numerical 




(21) 




(22) 
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Figure 4: The universal scaling plots of the order parameter and its fluctuations 
(inset) for D = 5 and for various values of the external field {h = 510~^, 710~^, 
10~^, 710"''). The dashed lines correspond to the mean-field solutions RMpix,!) 
and Dmf{x, 1) [see Eqs. (^^1^ ]. 

data and R{x, 1) is obtained. This is a confirmation of the RG-result D^ = A |441l45j . 
To the best of our knowledge no numerical evidence that five-dimensional DP 
exhibits mean-field scaling behavior was published so far. 

The rescaled fluctuation data is presented in Fig. HI As for the universal order 
parameter, the data of the fluctuations are in agreement with the corresponding 
universal mean-field scaling function 

L»mf(x,1) = 1 + . (23) 

2 V'l + {x/2f 

4.3 At the upper critical dimension 

At the upper critical dimension = 4 the scaling behavior is governed by the 
mean-field exponents modified by logarithmic corrections. For instance the order 
parameter obeys in leading order 

p^{dp,h = 0) oc dp\lndp\^, (24) 
p,{5p = 0,h) oc Vh\lnh\^ . (25) 

The logarithmic correction exponents B and T, are characteristic features of the 
whole universality class similar to the usual critical exponents. Numerous theo- 
retical, numerical, as well as experimental investigations of critical systems at 
have been performed (see for instance [MlllZlllHlllilSnilinillZlIlHlIllElEl)- 
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Logarithmic corrections make the data analysis quite difficult. Hence most investi- 
gations are focused on the determination of the correction exponents only, lacking 
the determination of the scaling functions at . 

Recently, a method of analysis was developed to determine the universal scaling 
functions at the upper critical dimension [^. In this work the authors use the 
phenomenological scaling ansatz (all terms in leading order) 

a,p.{6p,h) ~ A-^MF|lnA|'^(ap5pA|lnA|^a;,/lA'"MF|lnA^, (26) 

with Pmf = 1 and (Tmf = 2. Therefore, the order parameter at zero field {h = 0) 
and at the critical density {6p = 0) are given in leading order by 

a,p.{Sp,h = 0) ~ a^6p\lnaM^ R{1,0), (27) 

a,p.{5p = 0,h) ~ v^llnv^l^ ^(0,1) (28) 

with B = 6 + / and S = s/2 + I. Similar to the case D ^ the normalization 
i?(0, 1) = i?(l,0) = 1 was used. According to the ansatz Eg. (|26]) the scaling 
behavior of the equation of state is given in leading order by 

a^p^iSp,h) ~ yaji I In -R(x, 1) (29) 

where the scaling argument is given by 

X = ap5p^/aji I In ^ a^h |" (30) 

with H = 6 - s/2 = B - S. 

In case of directed percolation it is possible to confirm the scaling ansatz Eq. H26|) 
by a RG-approach j23j . In particular the logarithmic correction exponents are given 
by I = 7/12, h = —1/4, and s = —1/2. Thus the scaling behavior of the equation 
of state is determined by the logarithmic correction exponents [231 

B = S = 1/3, H = 0. (31) 

It is worth mentioning that in contrast to the RG results below the upper criti- 
cal dimension the logarithmic correction exponents do not rely on approximation 
schemes like e- or 1/n-expansions. Within the RG theory they are exact results. 

Similarly to the order parameter the following form is used for its fluctua- 
tions pu] 

a^/\p.X5p,h) ~ A^'llnAl'^ L'(ap(5p A | In A|^ a^/i A"'" | In A|") . (32) 

Using the mean-field value 7' = and taking into account that numerical simula- 
tions show that the fluctuations remain finite at the critical point (i.e. A; = 0) the 
scaling function 

a^^pS5p,h) ~ D(x,l) (33) 
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Figure 5: The universal scaling plots of the order parameter and its fluctuations 
(upper left inset) at the upper critical dimension = 4: for various values of 
the external field {h = 510'^, 210'^, 810"^, 410"^, 210-^). The logarithmic 
correction exponents are given byB = S = l/3 [2SI and H = 0. The right insets 
show the order parameter at the critical density and for zero field, respectively. The 
order parameter is rescaled according to Eqs. H27I28() . Approaching the transition 
point {h ^ and 6p — > 0) the data tend to the function f{x) = x (dashed lines) 
as required. 

is obtained, where the scaling argument x is given by Eq. with H = 0. The 
non-universal metric factor is determined again by the condition D{0, 1) = 1. 

Thus the scaling behavior of the order parameter and its fluctuations at 
is determined by two exponents (B = 1/3 and S = 1/3) and four unknown non- 
universal metric factors (a^, a^, a^, Oa). Following [20] we determine these values in 
our analysis by several conditions which are applied simultaneously: first, both the 
rescaled equation of state and the rescaled order parameter fluctuations have to 
collapse to the universal functions R{x, 1) and D{x, 1). Second, the order param- 
eter behavior at zero field and at the critical density are asymptotically given by 
the simple function /(x) = x when plotting [a^p^{5p,0)/apSp]^/^ vs. | lnap6p\ and 
[«aPa(0, h)/^/aJ^]^/'^ vs. | In ^/aJ^\, respectively. Applying this analysis we observe 
convincing results for B = S = l/3, H = 0, and for the values of the non-universal 
metric factors listed in Tabled The corresponding plots are presented in Fig.|SJ 

5 UNIVERSAL AMPLITUDE COMBINATIONS 



In the following we consider several universal amplitude combinations (see |lj for 
an excellent review). As pointed out in these amplitude combinations are very 
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Figure 6: The universal scaling function C{l,x)/C{—l,x) for various dimensions. 
The dashed lines correspond to an e-expansion of a RG approach [221 • The universal 
amplitude C(l, 0)/C(— 1, 0) is obtained from the extrapolation a^hlap dpl'"^ 0. 



useful in order to identify the universality class of a phase transition since the am- 
plitude combinations vary more widely than the corresponding critical exponents. 
Furthermore, the measurement of amplitude combinations in experiments or sim- 
ulations yields a reliable test for theoretical predictions. In particular, estimates 
of amplitude combinations are provided by RG approximation schemes like e- or 
1 / n-expansions. 

Usually numerical investigations focus on amplitude combinations arising from 
finite-size scaling analysis. A well known example is the value of Binder's fourth 
order cumulant at criticality (see e.g. [HB). Instead of those finite-size properties we 
continue to focus our attention to bulk critical behavior since bulk amplitude com- 
binations are of great experimental interest. Furthermore, they can be compared 
to RG-results |221- 

The susceptibility diverges as 

X{6p>0,h = 0) ~ a^,+ 6p-\ (34) 
X{6p < 0, /i = 0) ~ a^,_ i-6p)-^ , (35) 

if the critical point is approached from above and below, respectively. The ampli- 
tude ratio 

xjSp >0,h = 0) ^ Q;,,+ 

is a universal quantity similar to the critical exponents, i.e., all systems belonging 
to a given universality class are characterized by the same value a^,+/a^__. This 
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can be seen from Eq. ©. Setting ap|5p|A = 1 yields 

xjSp > 0, h) ^ C{+l,x) 
Xi6p<0,h) C{-l,x) 



(37) 



with X = ai^h\ap5p\~'^ . Obviously this is a universal quantity for all values of 
the scaling variable x. In particular it equals the ratio a^^^/a^^_ for x ^ 0, i.e., 
vanishing external field. In general, universal amplitude combinations are related 
to particular values of the universal scaling functions. 

In Fig. El the universal susceptibility ratio Eq. (|37)) is shown for various dimen- 
sions. The corresponding data saturates for x — > 0. Our estimates for the amplitude 
ratios ^(+1, 0)/C'(-l, 0) are 0.033 ± 0.004 for D = 1, 0.25 ± 0.01 for D = 2, as 
well as 0.65 it 0.03 for D = 3. In case of five-dimensional DP the amplitude ratio 
is constant, as predicted from the mean- field behavior 

Cmf(+1; x) _ ^ ^^g^ 

C'mf(~1) x) 

for all X. The behavior of the ratio C{+1, x)/C{—l, x) for D < reflects the 
crossover from mean-field to non mean-field behavior. Far away from the transition 
point, the critical fluctuations are suppressed and the behavior of the system is well 
described by the mean-fleld solution [Eq. ()38() ]. Approaching criticality the critical 
fluctuations increase and a crossover to the D-dimensional behavior takes place. 

In the already mentioned work [22|) Janssen et al. calculated the steady state 
scaling behavior of DP within a RG approach. In particular they obtained for the 
susceptibility amplitude ratio 



C(-1,0) 



V288 144 3 



e + 0{e 



2^ 



(39) 



leading to -0.2030 . . . for L> = 1, 0.2430 . . . for L> = 2, 0.6441 . . . for L> = 3. Except 
for the unphysical one-dimensional result these values agree well with our numerical 
estimates. 

Furthermore the parametric representation of the susceptibility was derived 
in |22j and it is straightforward to calculate the universal ratio Eq. (|37|) . The results 
are plotted for various dimensions in Fig. 1^1 It is instructive to compare these 
results with the numerical data since the theoretical curve reflects the accuracy of 
the RG estimations of all three quantities, the exponent, the scaling function, as 
well as the non- universal metric factors. All quantities are well approximated for 
the three-dimensional model. In the two-dimensional case we observe a horizontal 
shift between the numerical data and the RG-estimates. Thus the RG-approach 
yields good estimates for the exponents and the scaling function but the metric 
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factors are of significantly less quality. For D = 1 the e^-approximation does not 
provide appropriate estimates of the DP scaling behavior. Thus higher orders than 
0{e'^) are necessary to describe the scaling behavior of directed percolation in low 
dimensions. 

Analogous to the susceptibility the universal amplitude ratio of the fluctuations 
is given by 

Ap^{6p>0,h) D{+l,x) 



Ap^{6p<0,h) D{-l,x) 



(40) 



with X = ahh\ap6p\~'^ . In the case of absorbing phase transitions this ratio diverges 
for vanishing field. For Sp < the order parameter fluctuations are zero (absorbing 
state) for vanishing field whereas the fluctuations remain finite above the transition 
{5p > 0). Thus absorbing phase transitions are generally characterized by 

Si±h^ ^ oo. (41) 
D(-1,0) 

In Fig.Elwe plot the fluctuation ratio [Eq. (|4()j) ] as a function of the scaling variable 
ahh\ap6p\~'^ for various dimensions. We observe in all cases that the fluctuation 
ratios diverge for x — > 0. Only the one-dimensional system exhibits a particular 
behavior characterized by the minimum of the corresponding curve. The origin of 
this behavior is that for D = 2,3 the universal scaling function D{x, 1) exhibits a 
maximum for x > 0, whereas for D = 1 it is located at x < (see Fig. 8 in jl7j). 
In the five-dimensional model we observe a perfect agreement with the mean- 
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Figure 7: The universal scaling function D{l,x)/D{—l,x) for various dimensions. 
The dashed line corresponds to the mean-field scaling behavior. For D = 2 and 
D = 3 the mean-field curves are shifted the by the factors 9.56 and 2.12 to the left. 
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field behavior 

Dmf{+1,x) ^ 1 + yi + 4x _^ l + 2x 
DmA-'^,x) - 1 + VI + 4x ^0 2x ' 

Surprisingly, the two- and three-dimensional data are also well approximated by 
this formula provided that one performs a simple rescaling (x i-^ a^x) which results 
in Fig.[7|in a horizontal shift of the data. We suppose that this behavior could be 
explained by a RG-analysis of the fluctuations. 

Similar to the universal amplitude ratios of the susceptibility and the fluctua- 
tions other universal combinations can be defined. Well known from equilibrium 
phase transitions is the quantity (see e.g. [3]) 

= rd,B^^^, (43) 

which is also experimentally accessible, e.g. for magnetic systems. Here, T denotes 
the amplitude of the susceptibility x in ^^ro field (x ~ T 5T~'^) and B is the 
corresponding amplitude of the order parameter M (M ~ B6T^). The factor 
describes how the order parameter M depends on the conjugated field H at ST = 

In case of directed percolation these amplitudes correspond to the values B = 
a^^(l,0), r = ajo;, (7(1,0) as well as d, = a'^ RiO,!)'^ where 6 = a/ (3. The 
normalizations i?(l,0) = -R(0, 1) = 1 yield for the amplitude combination 

Rx = (44) 
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Figure 8: The universal scaling function C(l,x) for various dimensions. The 
dashed lines correspond to an e-expansion of an RG approach |22]. The universal 
amplitude R^ is obtained from the extrapolation ahh{ap5p)~" — > 0. 
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which is obviously a universal quantity. In Fig. [HI the scaling function C{l,x) is 
plotted as a function of x = a,j/i (a^ 5p)^^ for D = 1, 2, 3. The corresponding data 
saturates for x ^ 0. Our estimates are = 0.60±0.04 for D = l,R^ = 0.72±0.04 
for D = 2, and = 0.86 it 0.08 for D = 3. Note that the error-bars reflect only the 
data scattering in Fig.|Hl In contrast to the amplitude (7(1, 0)/C'(— 1, 0) the data 
of Ry- are affected by the uncertainties of the exponent 7 and the uncertainties of 
the metric factors a^, a^. These uncertainties increase the error-bars significantly. 
The two- and three-dimensional data agree quite well with the RG-results Ry. = 
0.7244 ... for L> = 2 and i?^ = 0.9112 ... for L> = 3 |22]. In the one-dimensional 
model the e^-expansion yields again an unphysical result (i?^ = —3.927. . . ). 

6 CONCLUSIONS 

We considered the universal steady state scaling behavior of directed percolation 
with an external field in D > 2 dimensions. Our data for D=5 coincide with the 
mean field solution, confirming that D^. = 4 is the upper critical dimension. At 
Z?c we presented for the first time a numerical scaling analysis of DP including 
logarithmic corrections. Our results agree well with those of a recently performed 
RG approach [SB]- Apart from the scaling functions we also considered amplitude 
ratios and combinations for the order parameter fluctuations and the susceptibil- 
ity. A comparison with RG results reveals that higher orders than 0{e^) are 
necessary to describe the scaling behavior in low dimensions. 
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